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Abstract 

We compare the absolute gain photometric calibration of the PlanckfiAFl and HerscheljSPWS. instruments on diffuse emission. The absolute 
calibration of HFI and SPIRE each relies on planet flux measurements and comparison with theoretical far-infrared emission models of planetary 
atmospheres. We measure the photometric cross calibration between the instruments at two overlapping bands, 545 GHz / 500 pm and 857 GHz / 
350 pm. The SPIRE maps used have been processed in the Herschel Interactive Processing Environment (Version 12) and the HFI data are from 
the 2015 Public Data Release 2. For our study we used 15 large fields observed with SPIRE, which cover a total of about 120 deg^. We have 
selected these fields carefully to provide high signal-to-noise ratio, avoid residual systematics in the SPIRE maps, and span a wide range of surface 
brightness. The HFI maps are bandpass-corrected to match the emission observed by the SPIRE bandpasses. The SPIRE maps are convolved 
to match the HEI beam and put on a common pixel grid. We measure the cross-calibration relative gain between the instruments using two 
methods in each field, pixel-to-pixel correlation and angular power spectrum measurements. The SPIRE/HFI relative gains are 1.047 (±0.0069) 
and 1.003 (± 0.0080) at 545 and 857 GHz, respectively, indicating very good agreement between the instruments. These relative gains deviate 
from unity by much less than the uncertainty of the absolute extended emission calibration, which is about 6.4 % and 9.5 % for HFI and SPIRE, 
respectively, but the deviations are comparable to the values 1.4 % and 5.5 % for HFI and SPIRE if the uncertainty from models of the common 
calibrator can be discounted. Of the 5.5 % uncertainty for SPIRE, 4 % arises from the uncertainty of the effective beam solid angle, which impacts 
the adopted SPIRE point source to extended source unit conversion factor, highlighting that as a focus for refinement. 

Key words. Methods: data analysis 


1. Introduction 

The Planck^ and Herschef satellites have an interconnected his¬ 
tory, including their launch together in 2009. Although operating 
with different scientific goals across many frequencies, Planck 
and Herschel have in common two very similar passbands, 
545 GHz / 500 pm and 857 GHz / 350 pm. This redundancy is 
very important for complementary studies. For example, the 
combination of Planck large-scale and Herschel small-scale ob¬ 
servations is valuable to the study of the cosmic infrared back¬ 
ground anisotropies (e.g., Planck Collaboration XVIII 2011; 
Viero et al. 2013) and to understanding the interplay between in¬ 
terstellar medium structures and their environment (e.g., Planck 
Collaboration XXII 2011; Juvela et al. 2011). With the 2015 
Public Data Release 2 (PR2) of the Planck data and the grow¬ 
ing public availability of processed Herschel observations, this 
is an opportune time to address the important question of the 
compatibility of measurements carried out by both instruments. 


* Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci¬ 
entific consortia funded by ESA member states and led by Principal 
Investigators from France and Italy, telescope reflectors provided 
through a collaboration between ESA and a scientific consortium led 
and funded by Denmark, and additional contributions from NASA 
(USA). 

^ Herschel (http://www.esa.int/Herschel) is an ESA space 
observatory with science instruments provided by European-led 
Principal Investigator consortia and with important participation from 
NASA. 


The Planck High Frequency Instrument (HFI, Planck 
Collaboration VI 2014) is composed of a set of 52 bolometers 
observing the sky at six frequencies (100, 143, 217, 353, 545, 
and 857 GHz). It provided full-sky maps with an angular reso¬ 
lution ranging from 9.'7 to 4.'6. The Spectral and Photometric 
Imaging Receiver (SPIRE, Griffin et al. 2010)^ provided photo¬ 
metric capabilities at 250, 350, and 500/rm (PSW, PMW, and 
PLW bands, respectively) to the Herschel Space Observatory 
(Pilbratt et al. 2010). The three arrays contain 139, 88, and 43 
detectors with angular resolution of 18!'2, 24!'9, and 36/3, re¬ 
spectively, and a common field of view of 4' x 8'. Larger fields 
were scanned and maps made from the time-stream data us¬ 
ing software available in Version 12 of the Herschel Interactive 
Processing Environment (HIPE)"^. 

In this paper we compare the absolute photometric calibra¬ 
tion of the HFI and SPIRE instruments in the overlapping band- 
passes (Fig. 1), focusing on the relative calibration (the cross¬ 
calibration relative gain) for emission that is extended on the sky, 
i.e., the so-called diffuse emission. The relative offset calibration 
will be presented separately (Schulz et al., in preparation)^. 

The paper is organized as follows. In Sect. 2 and Appendix A 
we briefly summarize the calibration schemes of PlanckfiWi and 
HerscheljSPlPfE frequency maps with relevant details concern¬ 
ing beams, bandpasses, and absolute calibrators. In Sect. 3 we 
describe the fields that we selected to study the cross calibra- 


The SPIRE Handbook is available at 
http://herschel.esac.esa.int/Docs/SPIRE/spire_handbook.pdf. 
http: //herschel. esac. esa. int/hipe/ 

^ For a preview, see http://herschel. esac. esa. int/TheUniverseExpl 
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tion of the diffuse emission. In Sect. 4 we discuss the factors that 
have to be taken into account to compare the intensity of a given 
source of emission in the different bandpasses of the two instru¬ 
ments. We then carry out two independent assessments of the 
relative gain using maps of diffuse emission; a direct correlation 
analysis in Sect. 5 and a power spectrum comparison in Sect. 6. 
We conclude in Sect. 7. 


2. Absolute photometric calibration 

Following the tradition of infrared and submillimetre experi¬ 
ments, pipeline processed Planck and Herschel measurements 
are calculated, calibrated, and reported as surface brightnesses 
(or monochromatic flux densities) denoted adopting as a 
reference spectrum a power law of index -1 (i.e., v S^efiv) - 
constant). Thus 

1 jN(y)S(y)dy 

VO jN(y)/vdy ’ ^ ^ 

where vo is the adopted nominal frequency of the hlter. For the 
standard HFI and SPIRE pipelines, vo is chosen to be equal to 
545 and 857 GHz, and to correspond to wavelengths of 350 and 
500 fim, respectively. The parameter AZ(v) is the net spectral re¬ 
sponse in terms of absorbed power including the aperture effi¬ 
ciency, the hlter spectral response, and for surface brightness the 
spectral dependence of the beam across the band (see, e.g., de¬ 
tails in Griffin et al. 2013 and the SPIRE Handbook, Chapter 5). 
Its normalization is unimportant because the response is used 
only in ratios. Absolute calibration therefore requires precise 
knowledge of the instrument properties, particularly bandpasses 
and beams, and there are related uncertainties summarized here. 
See Appendix A for more details. 

Recovering the true monochromatic brightness (or hux den¬ 
sity), 5y„, from5y„ requires an instrument-dependent colour cor¬ 
rection that depends on the shape of the SED, which varies over 
the sky. Because we are only interested in the relative calibration 
of the two instruments, these colour corrections are not needed 
explicitly. However, because of the different shapes of the spec¬ 
tral response functions for the pairs of bandpasses of the two 
instruments we are comparing (Fig. 1) and the different values 
of Vo, a bandpass correction dependent on the shape of the SED 
is needed for the cross-calibration analysis (see Sect. 4). 

The Planck/HFl calibration scheme for the 2015 PR2 data 
is described in Planck Collaboration VIII (2016). The 545 and 
857 GHz channels were calibrated by comparing measurements 
of Uranus and Neptune flux densities with predictions based 
on models of their emissivities (the so-called Uranus ESA 2 
and Neptune ESA 3 models produced by Moreno 2010 updating 
Moreno 1998). The statistical uncertainty of the measurements is 
1.1% and 1.4 % at 545 and 857 GHz, respectively, and the model 
predictions have an absolute uncertainty of 5 %. Combining the 
statistical and systematic uncertainties linearly, the overall HFI 
absolute calibration uncertainty for PR2 data is 6.1 % and 6.4% 
at 545 and 857 GHz, respectively. 

The relative planet model uncertainty between these two HFI 
bands, however, is expected to be of order 2%. Combining the 


® The spectral energy distribution (the SED) is often denoted P in 
the case of surface brightness and Fy in the case of point source flux 
densities. In this paper we will use Sy for both, the distinction being 
clear from the context. The pipeline product that we denote 5v„ is called 
5 pin in the SPIRE Handbook. 
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Figure 1. Spectral response for the P/anck/HFI 545 and 857 GHz 
Alters (red) and the //ericAeZ/SPIRE PLW and PMW Alters 
(blue), respectively, each normalized relative to a maximum at 
unity. For SPIRE these include the spectral dependence of the 
beam across the band as is appropriate for extended emission 
(compare Eig. 5.16 to Fig. 5.5 in the SPIRE Handbook v2.5). No 
correction is needed for HFI. 


statistical uncertainties in quadrature (thus 1.8 %) with this rel¬ 
ative planet model uncertainty linearly, the relative uncertainty 
between the two bands would be about 3.8 %. 

As small bolometer to bolometer differences and varia¬ 
tions over the sky have a negligible effect on our results 
(Appendix A. 1), we used a Gaussian beam with a full width half 
maximum (FWHM) of 4.'83 and 4.'64 at 545 and 857 GHz, re¬ 
spectively (Planck Collaboration VII 2016). 

The spectral responses shown in Fig. 1 are part of the Planck 
data release. We do not consider their uncertainties in this study; 
reported uncertainties on planet colour correction factors (Mars, 
Jupiter, Saturn, Uranus, Neptune) derived from these spectral re¬ 
sponses are about 0.01 % and 0.014 % at 545 and 857 GHz, re¬ 
spectively (Planck Collaboration IX 2014). 

Planck Collaboration XXX (2016) found that, with respect 
to the very accurate planet-independent calibration of the 100 
and 143 GHz channels, the relative photomeAic calibration of 
the 545 GHz channel is -f 2.3 % + 1.6 % using the solar dipole 
and -1-1.5%+ 1.8% using the Arst two CMB acoustic peaks 
(Appendix A. 1). This suggests that the uncertainty of the planet- 
based absolute calibration of the 545 GHz channel is not as great 
as the value of 6.1 % cited above, of which 5 % arose from the 
absolute uncertainty of the planet model. 

The photometric calibration of the //ersc/reZ/SPIRE instru¬ 
ment is described in detail in Bendo et al. (2013) and in the 
SPIRE Handbook. Calibration on the ESA 4 model of Neptune 
(Moreno 2012) introduces a 4% systematic uncertainty from 
model predictions. The statistical uncertainty on SPIRE photom¬ 
etry of the calibrators is about 1.5 % (cf. Bendo et al. 2013) in 
the 350/tm (PMW) and 500//m (PLW) bands. These two contri- 
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butions add linearly to produce a 5.5 % uncertainty in the point 
source calibration. 

The pipeline processing used to produce the SPIRE maps 
needed in this paper made use of a conversion from point 
source flux density [Jy beam '] to extended source surface 
brightness [MJysr^'] (Griffin et al. 2013), again for a refer¬ 
ence spectrum v 5ref(v) = constant. In the notation used in 
the SPIRE Handbook v2.5 the conversion factor may be writ¬ 
ten as /fptoE = (^ 4 E/^ 4 p)/fieff to highlight the inverse depen¬ 
dence on Qeff, which is the so-called effective beam solid an¬ 
gle^ for that bandpass accounting for the frequency dependence 
of the beam and calculated for the reference spectrum. The ra¬ 
tio Kh^IKa,? of monochromatic conversion factors defined in the 
Handbook, also calculated for this reference spectrum, is close 
to unity. Eollowing Table 5.2 of the SPIRE Handbook v2.5 as 
incorporated in the HIPE Version 12, the values of this ratio and 
Qeff that we used are 1.0015 and 0.9993 and 822.58 arcsec^ and 
1768.66 arcsec^, so that ^fptoE is 51.799 and 24.039 MJysr^' 
Jy"' for PMW and PEW, respectively. Underlying spectral re¬ 
sponses of SPIRE for extended emission are shown in Eig. 1 . 

The uncertainty of /fptoE is dominated by the 4 % uncertainty 
of the effective beam solid angles (Griffin et al. 2013). Adding 
this to the 5.5 % uncertainty of the point source calibration, the 
total uncertainty of the extended source brightness calibration 
amounts to 9.5 %. Ongoing efforts to improve the modelling of 
the radial beam prohle are described in Appendix A.2 and the 
implication of any change is quantihed in Sect. 7. 

Although Neptune is the reference calibrator in common 
with both instruments (Uranus is also used for HEI), two dif¬ 
ferent versions of the ESA model have been used. As described 
in Appendix A. 3, the net effect of using the ESA4 model for 
the HEI calibration, combined with the absolute photometric 
calibration derived from Uranus model, would be an increase 
of HEI brightness of 0.31 % and 0.16% at 545 and 857 GHz, 
respectively. Compared to the estimated uncertainties in cross¬ 
calibration relative gains found below in Sect. 5.2, this system¬ 
atic uncertainty is not a major concern. 

3. Field selection and map preparation 

Erom the Herschel Science Archive (HSA)*" we selected 15 
helds covering an appropriate range of average surface bright¬ 
ness and dynamics: two Science Demonstration Phase Hi-GAL 
helds at h = 0° (L59 centred at / = 59° and L30 at / = 30°), 
ten other helds in the Galactic plane from the Hi-GAL survey 
(Molinari et al. 2010), and three helds at higher latitude; these 
are Aquila (Konyves et al. 2010), Polaris (Miville-Deschenes 
et al. 2010), and Spider (in the North Celestial Pole Loop; e.g., 
Martin et al. 2015). Table 1 lists the held name, centre coordi¬ 
nates, size, and average surface brightness in the two PlancklHYl 
bands of interest. 

The SPIRE observations were carried out by mapping the 
held with a series of scans in one (the nominal) orientation and 
then a second time (with a different obsid) with a nearly orthogo¬ 
nal orientation (the cross-scan). Maps were made by processing 

^ This solid angle results from integration of the beam 
over 47r steradians, as distinct from the main beam solid an¬ 
gle for the main lobe, and hence corresponds to what else¬ 
where is called the antenna beam solid angle (see, e.g., 
http://ipnpr.jpl.nasa.gov/progress_report/42-64/64T.PDF) 

® The coefficient of determination provides a measure of the fraction 
of variance of the data explained by the model. 

http://archives.esac.esa.int/hsa/aio/doc/ 


simultaneously the data from scans in both the nominal and or¬ 
thogonal orientations. Eor each field, we used the Level 0 SPIRE 
timelines products from the HS A. The observations are then pro¬ 
cessed with HIPE to produce PMW and PEW maps with pixel 
sizes 10" and 14", respectively. Destriping was applied using 
the HerschellElPE destriper module. A constant offset was also 
removed that is equal to the median flux level of the map, be¬ 
cause SPIRE (like HEI) only measures relative fluxes, but this is 
of no consequence for the relative gain comparison. 

The HEI maps of the same fields were produced as cutouts of 
the 545 and 857 GHz whole-sky HEALPix" maps, corrected for 
zodiacal light, from the 2015 Planck data release (PR2; Planck 
Collaboration VIII 2016). We reprojected them on the SPIRE 
maps coordinates with a pixel size of 3' x 3'. 

Because the SPIRE beams are much smaller than those of 
HEI, we convolved the SPIRE maps with the HEI effective 
Gaussian beams (Sect. 2). We checked that the approximation 
of ignoring slight beam variations over the sky has a negligible 
impact on our analysis. Einally, we interpolated the convolved 
SPIRE maps to the HEI pixel grid. 

3.1. Estimated noise maps 

To estimate the high-frequency statistical noise for the HEI ob¬ 
servations we took advantage of the redundancy of Planck ob¬ 
servations during a stable pointing period (referred to as a ring; 
see Planck Collaboration VIII 2014). Among the products in the 
PR2 release there are two independent HEI maps produced by 
splitting each ring into two equal duration parts, the so-called 
half-ring maps (Planck Collaboration VIII2016, Appendix A. 1). 
The noise is assessed using the difference of cutouts from the 
two half-ring maps, following the method described in Miville- 
Deschenes et al. (2002) and Miville-Deschenes & Lagache 
(2005) in which the noise is assumed to be stationary. The his¬ 
togram of the brightnesses in these difference maps is very well 
described by a Gaussian in all the fields (see, e.g., the fit in Pig. 2, 
lower left, for the HEI Polaris field). The final estimated noise 
map (e.g.. Pig. 2, upper left) is weighted using the coverage. We 
note that this estimate is consistent with realizations based on 
the PR2 variance map (Planck Collaboration VIII 2016). 

The procedure for SPIRE was similar. We made separate 
maps from data for scans in each of the two orthogonal orien¬ 
tations (the two separate obsids) separately, creating a difference 
map that was then weighted by the coverage. In some SPIRE 
fields the histograms for the difference maps showed distribu¬ 
tions that are slightly broader than the central Gaussian core 
with low amplitude tails at high and low values of brightness 
difference. These are mostly due to 1 // noise (see also Sect. 6.1 
below) that is not perfectly removed when a map is produced 
from data in only a single scan orientation, as opposed to both 
(this is the fundamental reason for this observational strategy). 
Varying spatial coverage between the orthogonal sets of obser¬ 
vations, which is especially apparent along the map edges, also 
contributes to broadening the histogram. However these effects, 
which are already small at the SPIRE angular resolution, become 
completely negligible when the fields are brought to the HEI an¬ 
gular resolution. Similarly, the values shown here for the SPIRE 
difference maps are not directly comparable to the noise esti¬ 
mates discussed in Sect. 5.1 because the HEI 545 and 857 GHz 
beams encompass roughly 450 and 870 pixels in the correspond¬ 
ing SPIRE map, respectively. At the same resolution and pix- 


" See Gorski et al. (2005) and http://healpix. sf.net 
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Table 1. SPIRE fields selected for the SPIRE-HFI cross-calibration study. 


Name 

Centre (RA, Dec) 
[degrees] 

Size 

[degrees] 

5v(545 GHz) 
[MJysr-*] 

5v(857GHz) 

[MJysr-*] 

Hi-GALL59. 

(295.437, 23.0302) 

2.5 X 2.5 

28.31 

86.86 

Hi-GALL30. 

(281.541, -2.6095) 

2.5 X 2.5 

63.53 

219.03 

Hi-GAL Field 0 0_ 

(266.422, -28.9363) 

2.5 X 2.5 

73.56 

264.37 

Hi-GAL Field 2_0_ 

(267.710, -27.0518) 

2.5 X 2.5 

75.65 

284.66 

Hi-GAL Field 4 0_ 

(268.921, -25.1574) 

2.5 X 2.5 

74.53 

270.44 

Hi-GAL Field 6_0_ 

(270.130, -23.2524) 

2.5 X 2.5 

66.18 

237.39 

Hi-GAL Field 8 0_ 

(271.302, -21.3456) 

2.5 X 2.5 

66.49 

236.40 

Hi-GAL Field 33_0 . . . 

(282.893, 0.0483) 

2.5 X 2.5 

95.67 

323.82 

Hi-GAL Field 35 0 . . . 

(283.896, 2.0066) 

2.5 X 2.5 

109.95 

360.11 

Hi-GAL Field 39_0 . . . 

(285.909, 5.9450) 

2.5 X 2.5 

95.41 

306.30 

Hi-GAL Field 332 0 . . 

(243.747, -50.9331) 

2.5 X 2.5 

144.86 

508.12 

Hi-GAL Field 334_0 . . 

(246.174, -49.3698) 

2.5 X 2.5 

130.04 

452.62 

Aquila . 

(277.427, -2.7803) 

3.7 X 3.9 

29.25 

87.14 

Polaris . 

( 51.862, 88.5574) 

3.5 X 3.5 

2.13 

6.14 

Spider . 

(159.490, 72.9949) 

4.5 X 4.2 

0.65 

1.81 

GHz - Half-ring map 

difference 


SPIRE PLW 

obsid differe 


Polaris 


I—^-1 
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I--1 
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Figure 2. Left: HFI 545 GHz noise map of the Polaris field obtained from the difference of the two independent half-ring maps 
of the same field (see text). The bottom panel shows the normalized histogram of brightnesses in the difference map (black) and 
a Gaussian fit to the distribution (red), with the mean p, dispersion cr, and coefficient of determination as an estimator of the 
“goodness of fit” Right: SPIRE PEW (500/im) Polaris field noise map and histogram obtained from the difference between maps 
produced from scans in the nominal and orthogonal orientation, appropriately reweighted by the coverage map. The individual 
3° X 3° and 3?5 x 3?5 subfields were mapped at different times and so the effective noise is lower where they overlap. Likewise, the 
largest noise occurs at the edges of the subfields where the coverage by the bolometer arrays is least. 
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SPIRE PEW obsid difference 



AS^ [MJy.sr"'] 

Figures. Same as Fig. 2 but for a SPIRE difference map dom¬ 
inated by systematic effects (Hi-GAL field 37_0, 2?5 square). 
Although these fields might potentially be recovered using 
specifically tuned pipelines, we decided to discard them for rea¬ 
sons of consistency and accuracy. 


elization the distribution of SPIRE noise has a much smaller 
variance than is found for the noise in the HEI map. 

The SPIRE Polaris held (Pig. 2, right, and Pig. 8 below) is 
a special case, since it is the result of mapping two overlapping 
subhelds'^ (hence two pairs obsids), and the effect of weighting 
by the coverage is clearly seen in the estimated noise map. 

Pinally, we note that several additional helds were initially 
part of the analysis but were discarded because they displayed 
strong residual effects. An example of such a held is shown in 
Pig. 3. As a result of these residual systematic effects, adding 
such helds would not improve the accuracy of the joint estimate 
of the cross-calibration relative gain. Similarly, lower brightness 
helds are not useful because the instrumental noise limits the 
accuracy of the comparison. 


4. Bandpass corrections 

To compare Hux density or brightness measurements from two 
different instruments, it is necessary to apply a correction that 
takes the differences in their respective net spectral responses 
into account. This bandpass correction should not to be confused 
with the standard colour correction that, for a given bandpass, 
converts the brightness obtained at the nominal frequency from 
a given SED to that from another SED. A bandpass correction 

The structural and statistical properties of the subheld extending to 
the upper left were studied by Miville-Deschenes et al. (2010) using a 
power spectrum analysis. 




Figure 4. Correction factors from measurements with HEI to 
what would be observed by SPIRE, A'hfi->spire (Eq. 2), for 
MBB spectra (Eq. 3) with a range of temperatures, Tbbi and 
emissivity indices, (3, at the two HEI frequencies. We decreased 
the resolution along the temperature axis by a factor of 10 com¬ 
pared to the tabulated values to a step of 1 K. The anti-correlation 
between Tbb and (3 can produce the same MBB SED shape 
across the bandpass and, hence, the same correction. 


converts, for a given SED, the brightness obtained at the nominal 
frequency from a given bandpass to the brightness that would be 
obtained from another bandpass. Because the SED varies both 
within a given field and from field to field, bandpass corrections 
are essential in preparing the data for assessment of the cross¬ 
calibration relative gain; a corollary is that there is no unique 
bandpass conversion factor between the uncorrected HEI and 
SPIRE data. 

Although the two pairs of HEI and SPIRE bandpasses are 
similar, they are nevertheless sufficiently different that bandpass 
corrections are required, especially in the case of the 545 GHz 
and 500 jim bandpasses which overlap by only about two-thirds 
(see Eig. 1). Only after the bandpass correction has been ap¬ 
plied can one compare the two independent brightness measure¬ 
ments to investigate whether the absolute calibrations are con¬ 
sistent within their uncertainties. Technically, this comparison 
gives the desired measurement of the cross-calibration relative 
gain, which is expected to be independent of the SED and held 
on the sky. 
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spectrum v 5ref(v) = constant (see Eq. 1). We need a bandpass 
correction ^rHFi->spiRE to convert the observed brightness from 
HFI (at vq^i ) to that which would be observed in the paired SPIRE 
bandpass (at vqs), so that the bandpass-corrected HEI map can be 
compared to the SPIRE map. Thus 




VOh 


>SPIRE 




vos ^N\{{v)Sydv" 


( 2 ) 


where A/h and Ns are the net spectral responses for extended 
emission for HEI and SPIRE, respectively (as shown in Eig. 1). 


Figures. Planck HEI 857 to 545GHz brightness ratio versus 
MBB spectrum temperature, Tbb, for three different values of 
the emissivity index (J3 = 1.6, 1.8, and 2.0). For our application, 
dust temperatures would be in the low end of the range, allowing 
some discrimination of T bb- 




1.190 


1.185 


1.180 « 
A 
I 

1.175 X 

1.1 70 

■ 1.165 


Figure 6. Top: HFI 545 GHz spanning the Polaris held. Bottom: 
^HEi->spiRE bandpass correction, ^^545 at 545 GHz (Eq. 4), for 
the same held. 


4.1. Definition 

As described in Sect. 2, Herschel/SPIRE and Planck/HFl each 
provides maps of monochromatic surface brightnesses 5y„, at 
their individual reference frequency vo, assuming a reference 


4.2. Evaluating the bandpass correction 

Equation (2) shows that ^fHFi->spiRE depends on the shape of the 
source SED, Sy. In rare cases Syis well known beforehand (e.g., 
for a calibration source) but in general its shape has to be esti¬ 
mated (or assumed). 

A simple and robust model of the SED of thermal emission 
by Galactic dust in the submillimetre range is a modihed black- 
body (MBB) spectrum with a temperature, Tbb, and emissivity 
index, y 8 , in the form 


S y OC By(TBB) X ■ (3) 

Typical Tbb and (3 values for diffuse dust emission, which is the 
dominant component in the selected maps, range from 10 to 40 K 
and 1.4 to 2.5, respectively (e.g. Planck Collaboration XI 2014; 
Paradis et al. 2010). 

We used Eqs. (2) and (3) to produce tables of the Xhfi->spire 
bandpass-correction factors 

„ S PLW , „ S BMW ... 

K 545 - - and ^857 - p- (4) 

5 545 GHz 5 857 GHz 

for a range of Tbb from 10 to 40 K with a step of 0.1 K, and a 
range of /? from 1.2 to 2.2 with a step of 0.05. These are shown 
in Fig. 4. 

Several approaches can be used to estimate Tbb and (3 from 
the data, one of which is to hx the emissivity index at an appro¬ 
priate value deduced from previous studies and deduce the tem¬ 
perature from the measured ratio Planck HFI 857 to 545 GHz 
brightness ratio (Fig. 5). This ratio could be obtained from stan¬ 
dard linear regression between the two entire maps, but this 
would not take local temperature variations into account, and 
these variations can be important. Alternatively, the temperature 
can be deduced from the brightness ratio at each pixel. This 
is implemented in the HI PE pipeline developed to determine 
SPIRE maps offsets, with 8 '-resolution Planck maps. 

We implemented a second option, exploiting full-sky maps 
of Tbb and [3 produced from multi-frequency data as part of 
the Planck products (Planck Collaboration XI 2014; Planck 
Collaboration ES 2013). Neither Tbb and (3 are individually im¬ 
portant as long as, together in Eq. (3), they describe the slope of 
the SED that are fit in the data. Thus at every HFI pixel we can 
readily estimate the bandpass correction by interpolation in the 
above tables of K, producing bandpass-correction maps. An ex¬ 
ample is shown in Fig. 6 for the Polaris field at 545 GHz. We note 
that the bandpass corrections are different from pixel to pixel, 
with a range that is large compared to the uncertainty that we 
hnd for the cross-calibration relative gain below. Thus accurately 
computing and applying bandpass corrections is essential. 
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5. Pixel-to-pixel comparison 

After bandpass correction with ^fHFi->SPiRE an observation of a 
field with HFI should yield brightnesses equal to those measured 
with SPIRE (but the latter is missing an offset). This appears 
to be the case, as illustrated in Fig. 7. Any difference in slope 
compared to unity should be consistent with (or lower than) the 
combination of reported absolute calibration errors. 


5. 1. Degree of correlation 

The Pearson correlation coefficient between the pixel bright¬ 
nesses in the two maps (5J, in map i, in MJysr^') is the ratio 
of their covariance to the product of their standard deviations cr, 
(not to be confused with statistical uncertainties), i.e. 


Phh-spire “ 


Cov(5”*, 5™) 

^HFl ^SPIRE 


(5) 


which is unaffected by potential difference between the two cal¬ 
ibrations. Therefore, we use the Pearson correlation coefficient 
to assess the degree of linear correlation between the SPIRE and 
HFI surface brightnesses in a map. In the presence of noisy vari¬ 
ables, the computed correlation coefficient, Phh-spire, always un¬ 
derestimates the true value (see Sect. 5.3). 

Estimated correlation coefficients for each field at the two 
frequency bands are reported in Table 2. They average to 0.9986 
(+0.08%) and 0.9985 (±0.07%) at 545 and 857 GHz, re¬ 
spectively. This high degree of linear correlation between the 
SPIRE and HFI signals, although expected, conhrms that the 
SPIRE / HFI relative gain does not depend on brightness and is 
a useful validation that a simple linear ht procedure can be used 
to derive the cross-calibration relative gain and the SPIRE map 
offsets (Sect. 5.2). 


5.2. Relative gain measurements 


We used scatter plots to estimate the relative gain (G^‘^‘“) be¬ 
tween the two instruments. Following the hypothesis of linearity 
as demonstrated in Sect. 5.1, the relative gain is the slope of the 
correlation, independent of the zero levels. However, the same 
procedure can be used to hnd and restore the offset that has been 
removed from the SPIRE map during the map-making. Thus, the 
relationship is 


S 


SPIRE 

V 


= X 5^ - O 


( 6 ) 


Because the noise in the HFI map is considerably larger than 
in the convolved SPIRE map (Sect. 3.1), in practice we carried 
out the unweighted ht to the data in the sense of a regression of 
“X on Y” rather than “Y on X” and then inverted the slope to 
correspond to as written in Eq. (6). 

An example of the SPIRE-HFI scatter plot and computed 
linear ht is shown in Fig. 7 for the L59 held of the Hi-GAL 
survey. This analysis has been carried out for the 15 selected 
helds for the two frequency bands. Results for the relative gains, 
the slopes of the linear hts, are given in Table 2. The aver¬ 
age and standard deviation of the estimated values of is 
1.047 (+0.0069) and 1.003 (+0.0080) at 545 and 857 GHz, re¬ 
spectively. 


maps, we can attribute the dispersion about the regression line 
to the HFI uncertainties alone. We can estimate the Icr uncer¬ 
tainty, from the dispersion of the residuals in the horizontal 
direction; 5™ - (5™" + 0)/G^“‘. For the example L59 held, 
these residuals are shown in the lower panels of Fig. 7. 

Under the assumption of Gaussianity of the noise in our 
maps, which is reasonable here as shown in Sect. 3.1, the 
Pearson correlation coefhcient allows us to make an alternative 
estimate, a posteriori, of the noise in the HFI maps. Assuming 
that the measured signal can be written as 5meas = ^tme + 5 noise, 
which is a reasonable assumption, the computed correlation co¬ 
efhcient reads 


-^HH-SPIRE 


1 + 


1 + 


- 1/2 


-'HFI-SPIRE 9 


(7) 


where Mi are Icr uncertainties of the brightnesses characterizing 
‘5noise- In our selected helds, the signal variance cr^pj ^ t+spiRE 
is much higher than that of the uncertainty, i.e. m^ cr? in 
both SPIRE and HFI maps. Furthermore, as discussed above, 
“spire “hfi- Thus Eq. (7) simplihes to 


= 1 + 


(T 

KV u 


- 1/2 


-^Hn-SPIRE • 


(8) 


Because M^pj c+hfi’ variance of the measured 

HFI brightness as a hrst order estimate of cr^pj (the variance 
of the true, noiseless signal). Furthermore, we can assume that 
Phh-spire = 1 on the RHS. Using the computed Pearson correla¬ 
tion coefficients p in Table 2 on the LHS, we solve Eq. (8) for 
mhfi, designating this a posteriori estimate as M^pj. 

In the last two columns (for the two frequency bands) of 
Table 2, we tabulate the percentage difference of these two es¬ 
timates for each held: 100 KwJijpj - M^pj)/M^pj|. As expected, we 
hnd these to be in very good agreement (better than 0.4 %). The 
overall consistency further demonstrates that the gain does not 
depend on brightness and is consistent with the HFI brightness 
uncertainties dominating over those of SPIRE. 


5.4. Contributions to the error budget forG^'^^^ 

The relative gains in Table 2 show that at 857 GHz the SPIRE 
and HFI photometric calibrations agree to within 1 %. However, 
at 545 GHz there is a 4.7 % discrepancy that is statistically sig- 
nihcant (6.8cr). Each instrument has absolute calibration errors 
that were discussed in Sect. 2. As we showed in Appendix A.3, 
the Neptune models used for calibration of the two instruments 
are very close and can only explain at most a 0.31 % depar¬ 
ture from unity. Furthermore, the statistical measurement error 
from each instrument combined quadratically amount to an un¬ 
certainty of 1.86% and 2.06% in the relative gain at 545 and 
857 GHz. However, the SPIRE effective beam solid angle has 
an uncertainty of 4% (Griffin et al. 2013), acting as a system¬ 
atic effect on the SPIRE photometric calibration and thus on 
the gain estimates as well. This dominates the PlanckIHerschel 
cross-calibration error budget presented here and could account 
for the observed 4.7 % discrepancy at 545 GHz; the potential re¬ 
vision of the value of the solid angle discussed in Appendix A. 2 
would reduce the discrepancy to 2.7 %. 


5.3. Comparison of estimates of the uncertainty in the HFI 
maps 

Considering that the uncertainties in the convolved SPIRE maps 
are very small in comparison to the uncertainties in the HFI 


6. Cross calibration with power spectrum 
measurements 

Power spectra offer a means of probing the spatial structure of 
diffuse emission. In our context, we can use power spectra to 
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545 GHz [Planck/HFI] 



40 60 80 100 120 140 160 180 

HFI [MJy.sr"'] 


Figure 7. Correlation and pixel-to-pixel comparison between HFI and SPIRE observations of the Hi-GAL L59 field at 545 GHz. 
Top left: HFI and SPIRE maps of the L59 field. Top right: SPIRE-HEI pixel-to-pixel correlation. The linear fit from the unweighted 
“X on Y” regression is plotted in red, with the 1-to-l line in black for reference. The offset removed from the SPIRE map during 
the map-making has been obtained from the fit and added back. Bottom right: Horizontal residuals from the fit, as a function of 
brightness. Bottom left: Map of the horizontal residuals, with a colour range that is ten times finer (centred on 0), to enhance the 
detail. 


search for any scale dependence of the relative gain. This scale 
dependence is not expected, but can be introduced artificially by 
the data processing (any filtering process and the map-making). 
It is therefore of particular interest to check the level of agree¬ 
ment between power spectra of bandpass-corrected maps pro¬ 
duced by PlanckIWFl and those from i/eric/ieZ/SPIRE. 


However, the HFI and SPIRE spatial resolutions dif¬ 
fer greatly and that of HFI limits the largest spatial fre¬ 
quency (lowest spatial scale) for comparison to about k - 
0.3 arcmin *. As a result of map pixelization effects, we actu¬ 
ally limit our power spectrum comparison to k < 0.1 arcmin *. 
Furthermore, the SPIRE fields are limited in sky coverage, 
which limits the smallest spatial frequency to about k - 
0.007 arcmin ’ (proportionately lower for the three larger fields). 
At those scales, variations in the power spectra are dominated 
by cosmic variance but because cosmic variance is the same for 
the two maps it need not be considered in the comparison. This 
is not a large range in scale but thanks to the number of fields 


the power spectrum comparison still allows us to check for any 
systematic scale dependence of the gain. 


6.1. Implementation 

We implement the estimation of the power spectrum based on 
the methodology presented in Miville-Deschenes et al. (2002) 
via the IDE FFT routine. In overview, the power spectrum of an 
image Sx,y whose Eourier transform is Sk^^ky is computed from 
the map of the amplitude Ak^^k, defined by 

Aky,ky.^Sky,ky.Slky^\Sky,Kf- ( 9 ) 


The power spectrum P(k) dk is an angular average of Ak^^ky be¬ 


tween k and k + dk, where k 



To compare the power spectra of the dust signal we need to 
take the noise contamination and the effect of beam convolution 
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Table 2. Estimates of the relative gains, and Pearson correlation coefficients, Phfi-spire, at 545 and 857 GHz. 


Name 

/-’Scat 

‘J‘545 

^Scat 

^857 

P545 

P857 


Am**” 

% of M*” 

Hi-GALL59. 

1.036 

0.997 

0.9986 

0.9982 

0.14 

0.18 

Hi-GALL30. 

1.054 

1.002 

0.9989 

0.9993 

0.11 

0.07 

Hi-GAL Field 0 0_ 

1.058 

1.021 

0.9984 

0.9976 

0.16 

0.24 

Hi-GAL Field 2_0_ 

1.046 

1.007 

0.9992 

0.9992 

0.08 

0.08 

Hi-GAL Field 4 0_ 

1.039 

0.987 

0.9985 

0.9983 

0.15 

0.17 

Hi-GAL Field 6_0_ 

1.048 

1.006 

0.9988 

0.9990 

0.12 

0.10 

Hi-GAL Field 8_0_ 

1.047 

1.005 

0.9985 

0.9984 

0.15 

0.16 

Hi-GAL Field 33 0 . . . 

1.047 

1.000 

0.9993 

0.9991 

0.07 

0.09 

Hi-GAL Field 35_0 . . . 

1.045 

1.004 

0.9987 

0.9983 

0.14 

0.17 

Hi-GAL Field 39 0 . . . 

1.053 

1.005 

0.9995 

0.9994 

0.06 

0.06 

Hi-GAL Field 332_0 . . 

1.052 

1.006 

0.9981 

0.9970 

0.19 

0.30 

Hi-GAL Field 334 0 . . 

1.048 

1.005 

0.9989 

0.9991 

0.11 

0.09 

Aquila . 

1.047 

0.997 

0.9988 

0.9982 

0.13 

0.19 

Polaris . 

1.030 

0.989 

0.9983 

0.9981 

0.17 

0.19 

Spider . 

1.049 

1.010 

0.9958 

0.9978 

0.43 

0.22 

Average . 

1.047 (±0.0069) 

1.003 (±0.0080) 

0.9986 (±0.08%) 

0.9985 (± 0.07 %) 




Notes: 

Averages (and standard deviations) are over all fields, calculated a posteriori. 

Last two columns show the difference, in percentage, between the two HFI uncertainty estimates, as described in Sect. 5, which highlights that our 
assumptions about the FIFI uncertainties and the linearity between instrument measurements are very reasonable. 


into account. Formally, 


Psignal(^) ~ 


P{k) - Kik) 
B(k) 


( 10 ) 


where bn(k) is the noise power spectrum and B{k) the beam 
power spectrum. The noise power spectrum is computed from 
the estimated noise maps produced from redundant observations, 
as described in Sect. 3.1. As mentioned, there is a 1// noise 
component in the SPIRE estimate, which is clearly visible in 
the noise power spectra in Fig. 8. This component is not in the 
SPIRE power spectra computed on the total map, but because 
this component is about 20 - 50 times smaller than signal power 
spectrum at low k and it does not impact our analysis. 

In Fig. 8 we show how the deconvolved SPIRE power spec¬ 
trum at full resolution compares to that from HFI. This compar¬ 
ison is qualitatively interesting; as in the pixel-to-pixel compar¬ 
ison, we see a very good agreement across the range of spatial 
frequencies common to both instruments. Furthermore, the same 
power-law exponent appears to apply at the higher spatial fre¬ 
quencies accessible only to SPIRE, which is consistent with the 
study by Miville-Deschenes et al. (2010). 


operation in the image, which corresponds to a convolution of 
the true power spectrum in Fourier space, is carried out on both 
SPIRE and HFI maps, and so does not affect the power spectrum 
ratio and thus the relative gain determination. 

The relative gain is the square root of ratio of the two power 
spectra, i.e. 


G^^^ik) = 


Pspire(^) 

Phei(^) 


( 11 ) 


The relative gain is shown in Fig. 9 for each field, where it is 
seen that there is good agreement from field to field. The av¬ 
erage (and rms) over the spatial scales analysed is 1.043 
(+ 0.012) and 0.995 (+ 0.015) at 545 and 857 GHz, respectively. 
These estimates are in very good agreement (better than Icr) with 
the pixel-to-pixel average gains in Table 2. 

To search for any variation with spatial scale, we calculated 
the average G^^^(k) in five bins of k across the range 0.007 - 
0.1 arcmin *. For both frequencies we found the relative gains 
G^^^{k) to be stable with scale with rms variations lower than 
about 1.5 %. 


6.2. Relative gain from power spectra 

To compare SPIRE and HFI photometric calibrations quantita¬ 
tively, we compute P{k) on the bandpass-corrected HFI maps 
and on the SPIRE maps brought to the HFI resolution. This mit¬ 
igates against errors due to beam accuracy and differences in 
pixel window functions. 

The HIPE pipeline does not optimize the orientation of the 
final image to be used in the power spectrum analysis and so our 
maps feature zero coverage on the corners. Rather than rotat¬ 
ing the WCS, which would lead to significant distortion at small 
scales in the power spectrum, or cropping our data to small full- 
coverage insets (which was sufficient for the purposes of Fig. 8), 
which would reduce the range of overlapping scales, we zero- 
averaged our images and padded them with zeros. This masking 


7. Conclusion 

We performed a detailed comparison of the HerschellSPlRE 
and Planck/HFl absolute photometric calibration. We produced 
SPIRE maps for 15 fields that contain bright diffuse emission us¬ 
ing publicly released Herschel data. These maps were corrected 
for bandpass differences, reprocessed to match the HFI angular 
resolution and pixelization, and then compared to their Planck 
counterparts from publicly released data. The comparison was 
carried out using two methods: a pixel-to-pixel comparison and 
a power spectrum analysis. We discarded several fields that were 
originally in the study to limit systematic effects arising from 
significant SPIRE processing residuals and to maintain a high 
signal-to-noise ratio in the maps and a high accuracy in the com¬ 
parison. 


9 




















Bertincourt et al.: Planck[t{F\ and /feri-c/ieZ/SPIRE Cross Calibration 


POLARIS [SPIRE / PLW] 




0.001 0.010 0.100 1.000 10.000 
k [arcmin"'] 


POLARIS [SPIRE / PMW] 




0.001 0.010 0.100 1.000 10.000 
k [arcmln"'] 


Figure 8. Left: SPIRE PLW and PMW maps of the Polaris field. The black square outlines the region of the map that is used to 
compute the power spectra shown on the right. Right: Power spectra of SPIRE map (blue) and bandpass-corrected HEI (red) map 
for PLW (top) and PMW (bottom) in the Polaris field. The solid lines are the dust power spectra obtained from the image power 
spectra (dashed line), from which the noise spectra (dotted line) are subtracted followed by division by the effective beam power 
spectra (dash dotted line; arbitrarily shifted on the y-axis). 


The pixel-to-pixel comparison revealed the expected very 
high degree of linearity between the two datasets, and allowed 
for a robust estimate of the relative gain between the two instru¬ 
ments in each of the overlapping bands. The power spectrum 
analysis provided gain estimates with similar values and com¬ 
parable statistical uncertainties. The good agreement of the rel¬ 
ative calibration measured on distinct fields shows robustness to 
a number of potential systematic errors in either instrument, in¬ 
cluding nonlinear response, uncorrected time variability in the 
gain, and variation of the effective beam. 

We found that the estimated relative gains are well inside 
the absolute photometric uncertainties quoted for the two instru¬ 
ments of about 6.4% and 9.5% for HEI and SPIRE, respec¬ 
tively. However, because both instruments use the same plane¬ 
tary calibrator, that contribution to the relative uncertainty can 
be reduced to the very small difference between the ESA 3 and 
ESA 4 models of Neptune. Then the deviations are comparable 
to the remaining uncertainty values of 1.4% and 5.5 % for HEI 
and SPIRE. The difference in the ESA 3 and ESA 4 Neptune 
models used in the HEI and SPIRE calibrations, respectively. 


decreases the two relative gain estimates by only a very small 
fraction (0.31 % at 545 GHz and 0.16 % at 857 GHz). 

At 545 GHz the departure of the relative gain from unity is 
statistically signihcant, which raises the question of whether a 
systematic calibration correction should be applied prior to any 
comparative analysis between HEI maps at 545 GHz and SPIRE 
maps at 500jum. However, this departure could be accounted for 
by the 4 % systematic uncertainty of the SPIRE beam solid an¬ 
gle. Because this dominates the error budget of the SPIRE-HEI 
cross calibration, we provide here the SPIRE/HEI relative gains 
with their explicit dependence on the SPIRE unit conversion fac¬ 
tor /fptoE as follows: 


G 545 = 1.047 (+0.0069) x| 

1 -^fptoE \ 

124.039/ 

(12) 

Gs 57 = 1.003 (+0.0080) x| 

1 KptoE 

151.799/ 

(13) 


at 545 and 857 GHz (or PLW and PMW), respectively, where 
the normalizations of /fptoE are our adopted values in MJy sr ' 
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Figure 9. Relative gain, from the square root of the SPIRE-to-HFI power spectrum ratio, for each field at 545 GHz (left) and 
857 GHz (right). The mean gain from the power spectrum analysis (solid line) is in close agreement with the mean gain 
from pixel-to-pixel analysis (dotted line), which is well within the Icr dispersion of the latter (shaded in grey). 


Jy“'. These formulae would be useful for any future analysis that 
might make use of an updated version of the SPIRE ^fptoE (see 
e.g. Appendix A.2). 

We used two methods of comparison, both based on ex¬ 
tended emission, but did not make a comparison based on point 
sources. There are SPIRE counterparts for many extragalactic 
sources of the Planck point source catalogue (e.g., in the fields of 
H-ATLAS; Eales et al. 2010), but confusion noise at the Planck 
HFI channels induces uncertainties in point source photometry 
that are too large for an accurate estimate of the cross calibra¬ 
tion of the two instruments, especially compared to the level of 
precision that can be achieved with maps of diffuse emission. 
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Appendix A: Notes regarding absoiute photometric 
caiibration 

The absolute calibration of the two instruments is described in 
Sect. 2. Some supplementary details are provided here. 

A.1. Planck/HFI 

HFI is a scanning instrument and its effective beams relevant 
to the frequency maps are the convolution of (i) the optical re¬ 
sponse of the telescope and feeds; (ii) the processing of the time- 
ordered data and deconvolution of the time response; and (iii) the 
merging of several surveys to produce frequency maps. While 
accurate effective beams can be recovered at each sky position, 
this is a time-consuming process. However, the rms variations 
of the effective beam solid angle across the sky are 0.79 % and 
0.49 % at 545 and 857 GHz, respectively (see Table 3 in Planck 
Collaboration VII2016), which are negligible with respect to the 
calibration uncertainty, and so ignoring the spatial variation has 
a negligible effect for our study of diffuse emission calibration. 
Therefore, for convenience of implementation in our processing 
pipeline we decided to use an effective Gaussian beam instead, 
as cited in Sect. 2. 

Accurate knowledge of the spectral response is essential 
for comparing surface brightness obtained through different 
filters or for computing colour corrections. Spectral responses 
for PlanckfHFl were measured pre-flight for each detec¬ 
tor and checked with ground and in-flight data (see Planck 
Collaboration IX 2014). As noted in Sect. 2, these are known 
with sufficient accuracy not to impact the results in this paper. 

The photometric calibration for the PlanckfHFl channels 
which are CMB dominated (100, 143, and 217 GHz) or have 
a strong enough CMB signal (353 GHz) is achieved through ob¬ 
servations of the dipole arising from the orbital motion of the 
Planck spacecraft around the Sun (the orbital dipole). The cali¬ 
bration accuracy is very high, 0.09% and 0.07% at 100 and 143 
GHz, respectively (Planck Collaboration VIII2016). 

The 545 GHz channel and even more so the 857 GHz chan¬ 
nel cannot be calibrated on the orbital dipole, which is too 


weak with respect to Galactic dust emission, and so, as de¬ 
scribed above, these were calibrated using Uranus and Neptune. 
Nevertheless, given the high absolute calibration accuracy at 
100 and 143 GHz, the relative photometric calibration of the 
545 GHz channel with respect to those two channels can be used 
to assess the accuracy of the 545 GHz planet-based absolute cal¬ 
ibration. 

Planck Collaboration XXX (2016) have used the dipole gen¬ 
erated by the solar system motion with respect to the CMB (the 
solar dipole) and the first two acoustic peaks of the CMB fluctu¬ 
ations to measure the relative calibration of all Planck HFI chan¬ 
nels, except 857 GHz, adopting the average of 100 and 143 GHz 
as reference. As summarized in Sect. 2, it was found that the 
planet-calibrated data at 545 GHz has a relative calibration of 
-t-2.3 % + 1.6 % using the solar dipole and -1-1.5 % + 1.8 % using 
the first two CMB peaks. This is remarkable and suggests that 
the uncertainty of the planet-based absolute calibration of the 
545 GHz channel is not as great as the value of 6.1 % cited in 
Sect. 2. Recall that 5 % arose from the estimated absolute un¬ 
certainty of the planet model predictions, (Planck Collaboration 
VIII 2016) which relates mostly to uncertainties on the thermal 
profiles from Lindal (1992). Even if the uncertainty were that 
large, the relative planet model uncertainty between these two 
HFI bands is expected to be of order 2 %. 

The previous Planck 2013 data release (PRl) also relied 
on the Neptune and Uranus models for absolute calibration. 
Improvements in the data processing and calibration proce¬ 
dure shifted the final absolute calibration of frequency maps by 
-1.8 % and -3.3 % at 545 and 857 GHz, respectively (the 2015 
PR2 maps are fainter). 


The uncertainty of the conversion factor XptoE (Sect. 2) is domi¬ 
nated by the 4 % uncertainty of the effective beam solid angles, 
Qeft- We adopted values from Griffin et al. (2013) (HIPE Version 
12 ). 

Ongoing efforts to improve the modelling of the radial beam 
profile have revised the values of Qeff to 1804.3 + 13arcsec^ 
(+ 0.7 %) and 831.4 + 3.8 arcsec^ (± 0.5 %) for PLW and PMW, 
respectively. These amount to increases of just 2.0 % and 1.1%; 
furthermore, we note that the stated uncertainties in Oeft are now 
considerably less than the 4 % cited above. Adopting the revised 
values of Qeff would decrease the above conversion factors XptoE 
and thus decrease the brightness of the SPIRE maps and the 
SPIRE/HEI relative gains estimated in this paper (Eqs. 12 and 
13). 

Although these values have been incorporated in HIPE 
Version 14'^, nevertheless, until such time as revised beam solid 
angles are adopted officially in the SPIRE Handbook, and to fa¬ 
cilitate tracking of dependencies and changes, we felt it prudent 
to frame our results in terms of the cited values from the version 
of HIPE that we used (Version 12), consistent with the values 
from the latest SPIRE Handbook, v2.5. 


As seen in Eig. A.l, the main differences between the ESA 3 
model of Nepture used for the HEI calibration, and the ESA 4 
model of Nepture used for the SPIRE calibration rest in an up¬ 
dated treatment of the CO absorption features and HCN emis¬ 
sion lines. Several of these features fall within the HEI and 

http://herschel.esac.esa.int/twiki/bin/view/Public/SpirePh 
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Figure A.l. Top: Comparison of the ESA 3 (dashed black) and 
ESA 4 (green) Neptune models. Corresponding flux densities in 
the 545 and 857 GHz HEI bands are shown as triangles (ESA 3) 
and squares (ESA 4). Bottom: Relative difference of the two 
model spectra and the flux densities. 


SPIRE bands and so affect estimates of the cross-calibration fac¬ 
tor as a systematic error. We computed HEI equivalent Neptune 
flux densities for both models and find that, if calibrated on 
ESA 4 rather than ESA 3, the HEI absolute calibration factors 
on Neptune would increase by 0.63% and 0.31 % at 545 and 
857 GHz, respectively (see bottom panel of Eig. A.l). Combined 
with the absolute photometric calibration of HEI derived from 
Uranus model, these drop down to an increase of 0.31 % and 
0.16% at 545 and 857 GHz, respectively. We note that an in¬ 
crease of the HEI calibration factors translates into an increase 
of HEI brightness and thus a decrease of the SPIRE / HEI relative 
gains that we estimate in this paper (Eqs. 12 and 13). 
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